Emission properties of an oscillating point dipole from a gold Yagi-Uda nanoantenna array 
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We investigate numerically the interaction of an oscillating point dipole with a periodic array of optical Yagi- 
Uda nanoantennas in the weak coupling limit. A very strong near-field enhancement of the dipole emission by 
the resonant plasmon mode in the feed element is predicted in this structure. It is shown that the enhancement 
strength depends strongly on the dipole position, the direction of the dipole moment, and the oscillation fre- 
quency. The radiative intensity of the point dipole from appropriate places next to one feed element may exceed 
the radiative intensity of an equivalent dipole in free-space by a factor of hundred. In spite of only one director 
used in each nanoantenna of the array, the far-field emission pattern is highly directed. The radiative efficiency 
(the ratio of the radiative to the full emission) appears to be around 20%. 

PACS numbers: 78.67.-n, 42.60.Da, 73.20.Mf 



Nanophotonics has been in the focus of intensive inves- 
tigations in recent years. One of its numerous areas is the 
nanoscale controlling of light emission from a single molecule 
or quantum dot. Promisingtools for the realization of this goal 
are optical nanoantennas llll4l3l. 

Metal antennas are traditionally used for controlling the ra- 
diation pattern of electromagnetic wave emission in the radio 
and microwave frequency range. Though the electromagnetic 
properties of metals in the optical range differ significantly 
from that in the radio and microwave range, it seems to be 
reasonable to use the main concepts of radio antennas in the 
optical range as well. It has been suggested 10, [sl [sl |7, 13] 
to construct a nanooptical antenna with elements that are ar- 
ranged as in the radio Yagi-Uda antennas. 

Yagi-Uda antenna consists usually of one or two reflectors, 
one feed element and several directors with appropriately se- 
lected scattering phases (reflector and director are slightly de- 
tuned inductively and capacitively). As it has been recently 
shown, the nanoantenna elements can be nanorods llsils, 131. 
core-shell [2, 7] or spherical lll2ll nanoparticles. All elements 
scatter the light and the resulting interference forms a highly 
directed beam along the antenna axis. The size of the elements 
in such optical Yagi-Uda nanoantennas have to be smaller 
than the light emission wavelength in free-space. Such op- 
tical nanoantennas only work efficiently in narrow frequency 
domains, where the interaction of the emitter with light is res- 
onantly enhanced because of the excitation of so-called local- 
ized plasmons lfl4l[l5ll in the nanoantenna elements. 

Spontaneous emission is not an intrinsic atomic property, 
but it depends sensitively on the local density of photonic 
modes at certain frequencies in a microcavity iflH 171. or, 
equivalently, on the local electror nagnet ic field value at the 
position of the quantum emitter 13 U 111- Using resonances, 
it is possible to increase the local electromagnetic field sig- 
nificantly, and, resultantly, to enhance and redirect the dipole 
emission. In the case of localized plasmon resonances, the 
collective excitation of electrons at the plasmon frequency 



leads to a considerably enhanced emission rate when the point 
dipole is located in the vicinity of metallic nanoparticles with 
the appropriate orientation of its dipole moment Jsl S ■ 

The exact description of photon radiation from the quantum 
emitter located in some metal-dielectric environment is very 
complicated. A convenient approximation is a model of an 
oscillating point dipole. It oscillates with constant frequency 
and magnitude fixed by the external source (so-called weak 
coupling limit). In other words, the emission of a current 
j{r, t) - jo ■ 6{r - ro) • e^'"' inside an environment with spa- 
tially modulated permittivity has to be calculated. This system 
is now classical and can be described by Maxwell's equations. 

The goal of this paper is to calculate the radiation pattern 
and emission rate of one oscillating point dipole located in the 
periodic aiTay of optical Yagi-Uda nanoantennas. Each an- 
tenna in the array consists of three rectangular gold elements, 
director, feed, and reflector. The point emitter is located a few 
nanometers from the edge of one feed element. To the best 
of our knowledge, the emission properties of an oscillating 
point dipole from such a structure have not been investigated 
yet. The so far investigated structures have been either single 
Yagi-Uda nanoantennas S S Eltl, or arrays of simpler 
spherical shapes Using arrays of antennas in combina- 
tion with several emitters with controlled phase difference of 
oscillation, we can control the emission directivity addition- 
ally llSll . as in phased antenna arrays. 

The calculation of the radiation characteristics from oscil- 
lating dipoles in two-dimensional photonic crystal slabs can 
be performed using the method of scattering matrix |3 2oll . 
Its main idea is to divide the structure into two parts adjacent 
to a plane through the dipole, to find the scattering matrices of 
these parts, and to combine them by considering the bound- 
ary conditions at the plane with the dipole. We use the Fourier 
modal method (FMM) IH, 21 1, which can be improved by the 
formulation of the correct Fourier factorization rules ll22ll and 
adaptive spatial resolution jJst], as well as of matched coordi- 
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FIG. 1: (Color online) Lateral view of a periodic array of optical 
Yagi-Uda nanoantennas with an x-oriented oscillating point dipole 
(black point with arrow) located at a horizontal distance of 5 nm 
from the feed element of one nanoantenna. 



nates for more complex metallic shapes Il24tl . 

The structure of interest is shown schematically in Fig. [T] 
It consists of a glass superstrate (e = 2.4), a periodic array of 
optical Yagi-Uda nanoantennas ordered in a rectangular lat- 
tice, and a quartz substrate (e - 2.13). The periods along 
the X- and y-axis are equal to 450 nm and 300 nm, respec- 
tively. Each antenna consists of three rectangular gold paral- 
lelepipeds of 30 nm height and 100 nm width. The lengths 
of the top (director), middle (feed), and bottom (reflector) el- 
ements are 220, 250 and 300 nm, respectively. They are lo- 
cated in glass and the vertical distance between them is equal 
to 100 nm. We assumed the gold permittivity to be described 
by the Drude formula with 9016 meV plasma frequency and 
8 1 meV damping rate. 

Our first goal is to calculate the directional pattern as well 
as the emission spectra in the direction normal to the antenna 
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FIG. 2: (Color online) Calculated spectra of the emission in top (blue 
curve) and bottom (green curve) directions of the oscillating point 
dipole (directed along x) coupled to the gold Yagi-Uda nanoantenna 
array. The intensity is normalized to the maximum radiation inten- 
sity of an equivalent dipole (with the same magnitude and the same 
frequency) in free space. 



plane. It makes sense to normalize the computed emission 
intensity P(i?, tf), i. e., the Poynting vector of the dipole emis- 
sion in the far field as a function of the spherical angles and 
(f), to the maximum intensity of the emission of a point dipole 
in free-space, which oscillates with the same magnitude and 
frequency. Thus, we can easily distinguish the enhancement 
of emission (P > 1) from the attenuation (P < 1) compared to 
the dipole located in homogeneous vacuum. 

To characterize quantitatively the full dipole emission and 
its radiative part, it is convenient to define the Purcell factor 
and its radiative part as 
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respectively. Here, So and 2 denote spheres with an infinitely 
small and large radius, respectively, that surround the point 
dipole. P is the Poynting vector of the dipole emission from 
the antenna aiTay, is its counterpart for dipole emission in 
free space. The antenna absorption losses can be characterized 
by the non-radiative part of the Purcell factor Fp, 



Fp-F, 



rad 
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Now, we can also introduce the radiative efficiency i] indicat- 
ing the radiative part of dipole emission, 



Ff 



(3) 



In what follows the results of these quantities will be presented 
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employing 1633 spatial harmonics in the FMM 1121 . 

The calculated emission spectra in top (+z) and bottom (-z) 
directions of the x-directed oscillating point dipole located in- 
side the periodic array of gold Yagi-Uda nanoantennas are 
shown in Fig. 121 The dipole is placed on the horizontal sym- 
metry axis along .^-direction at a distance of 5 nm to the edge 
of the feed element (see in Fig.[T]i. Only one strong and narrow 
resonance occurs in the emission spectra in the top direction 
at the photon energy fioj - 820 meV {A - 1.5 //m). Its mag- 
nitude is about 250 and the FWHM is 107 meV. The emission 
to the bottom direction is significantly smaller. It has even a 
minimum at the plasmon resonance. 

Figure [3] depicts the calculated radiation pattern of emis- 
sion P{-&, (f) Cr at the resonant photon energy fia) - 820 meV 
(where is the radial unit vector of a spherical coordinate 
system centered in the point dipole position). The calculated 
radiative part of the Purcell factor according to the second 
Eq. ([U appears to be as high as 80 (see also Fig. |4] below), 
which indicates a very strong enhancement of the dipole emis- 
sion by the feed element. 

In spite of only one director employed in our Yagi-Uda an- 
tenna, the emission in top direction appears to be highly di- 
rectional. In order to characterize it, the angular directivity 
D{&, ip) can be calculated [|25tl , which indicates the part of the 
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FIG. 3: (Color online) Calculated 2D polar diagrams (solid lines) of far-field emission in xz (left panel) and yz (central panel) planes, as well 
as the full 3D-directional diagram (right panel) for the x-oriented oscillating point dipole coupled to the gold Yagi-Uda nanoantenna array at 
resonant photon energy {hoj = 820 meV). The emission intensity is normalized to the maximum radiation intensity of the equivalent dipole 
(with the same magnitude and the same frequency) in free space and shown in the right panel as colored surface P{&, tp) e,-. The dipole position 
is in the center-of-coordinates. Dashed lines in left and central panels show the scaled by a factor of 1/2 polar diagrams of P2(i,&,<p), the 
emission of a pair of synchronized dipoles separated by distance d (see explanation in the text). 



full emission radiated along the direction (p), 

4nP(-9, (fi) 
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Its maximum value Z)™^" = max[D(j?, ip)] is called directiv- 
ity and indicates the antenna's ability to form a narrow beam. 
The larger the directivity, the narrower is the light beam. An 
isotropic radiator would have a directivity of 1, whereas for 
an oscillating point dipole in free space, Z)™'' - 1.5. In the 
case of the radiation pattern of Fig. [3] Z)™'' - 4.7, i.e., it ex- 
ceeds the directivity of the dipole in free space by a factor of 
3.1. In jsl the directivity of a single Yagi-Uda antenna with 
three directors has been calculated as 6.4, which is only 1.37 
times larger than in our case with only a single director. Fur- 
thermore, we can also increase the directivity using more than 
one dipole emitter coupled to different antenna array elements 
which is shown below. 

It is also instructive to investigate the dependence of the 
emission enhancement on the dipole position and on the ori- 
entation of its dipole moment respectively to the feed ele- 
ment. Figure |4] shows the calculated dependence of the ra- 
diative part of the Purcell factor F^^'^ with respect to the dis- 
tance between the feed element and the dipole for x-, y-, and 
z-directed dipoles. The emission enhancement of the antenna 
array decreases with increasing distance. The strong polariza- 
tion and distance dependence demonstrates the local nature of 
the antenna-dipole enhancement. 

For the different distances indicated in the inset of Fig. |4] 
the emission is only significantly enhanced for the x-polarized 
dipole, which means that a plasmon mode with charges oscil- 
lating along X is excited in the system. Such mode leads to a 
very strong enhancement of emission and determines the di- 
rectional pattern. As the electromagnetic field near the edges 



of the feed element is known to have a dipolar character at 
the fundamental plasmon mode, it is not surprising that the 
radiation pattern of the z-polarized dipole displaced vertically 
by 10 nm above the edge of the feed element (see the inset 
in Fig.|5]l nearly coincides with that of the .x-polarized dipole 
5 nm apart from the edge of the feed element (see Fig. |3]l. 
Unlike a z-polarized dipole in free space, the majority of the 
emission is directed along the dipole polarization (i.e., along 
the z-axis). The magnitude, however, is about 3 times smaller 
than that of the x-oriented dipole. 

In this point let us discuss the results of the full Purcell 
factor Fp and the radiative efficiency rj. The calculated re- 
sults are shown in Fig. |5] where the dipole is placed 10 nm 
above the feed (i.e., in a homogeneous transparent glass). The 
Purcell factor is rather large and reaches the maximum of 
Fp X 200 - 300 when the dipole is located above the edge 
of the feed element. Our numerical analysis shows that the 
growing oscillations of the calculated Purcell factor when the 
,Y-coordinate approaches the edges of the director, feed, and 
reflector elements (vertical dashed lines in Fig. |5]i are the ar- 
tifacts of our calculation method. The radiative efficiency ap- 
pears to be almost independent of the dipole's x-coordinate 
and is about 20% when the emitter is above the feed element. 

It should be mentioned that we omitted the calculation of 
the full Purcell factor as well as the radiative efficiency 77 in 
Fig.m where the dipole is placed inside a modulated layer. In 
this case, it appears that we cannot calculate these quantities 
correctly in the FMM, because the discontinuous permittivity 
function of such a layer is described by a truncated Fourier 
expansion. Thus, due to the Gibbs phenomenon, also the non- 
absorptive surrounding of the metal exhibits a small imaginary 
part. The point dipole approximation, however, fails if the 
dipole is placed inside a lossy material. 
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Finally, we would like to mention that the antenna array 
allows to control the emission directivity further, by using 
several emitting dipoles in different positions and control- 
ling their phase difference of oscillations. Figure [3] shows 
the scaled by 1 /2 2D polar diagrams of the emission of two 
synchronized dipoles, positioned at neighboring antennas of 
the array at distance cf, proportional to integer numbers of 
periods (see dotted and dashed lines in the left and central 
panels). The emission intensity of the synchronized pair of 
dipoles normalized to the maximum emission intensity of 
two non-synchronized dipoles is then simply P2id,§,(p) - 
P('&,if)^l +cos(27m(cf •e"^)//l)j, where P{§,(p) is the direc- 
tional pattern of a single dipole (solid line), and n is the re- 
fractive index of superstrate. Note that the polar diagrams Po 
of the synchronized dipoles demonstrate the effect of super- 
radiance: the Purcell factor of two synchronized closely posi- 
tioned dipoles becomes twice larger than unsynchronized. It 
is seen also that it becomes possible to increase the directivity 
using two emitting dipoles. 

To conclude, we show that the gold Yagi-Uda nanoantenna 
array enhances and simultaneously directs radiation of the os- 
cillating point dipole. The dipole-nanoantenna enhancement 
depends very strongly on the oscillating frequency, dipole po- 
sition, and orientation of its moment. It becomes possible to 
control the emission directivity further by using several syn- 
chronized emitting dipoles attached to different antennas in 
the array. This opens a way to manipulate an excited-state 
lifetime of a quantum emitter and to fabricate narrow beam- 
ing nanoscale antennas in the optical range. 

We acknowledge support from BMBF (13N 9155, 13N 
10146), DFG (FOR 557, FOR 730, GI 269/11-1), DFH/UFA, 
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FIG. 4: (Color online) Calculated dependencies of the radiative part 
of the Purcell factor FjJ"' on the horizontal distance between the 
dipole and the edge of the feed element for a x-, y-, and z-polarized 
dipole (squares, circles and diamonds, respectively). The vertical 
dashed line marks the position exactly above the edge of the reflec- 
tor. The geometry is explained in the inset: the x-polarized dipole 
is shown as a black dot with arrow, it is is centered with respect to 
the feed element along y- and z-direction and shifted in x-direction, 
along the dotted horizontal line. 
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FIG. 5: (Color online) Calculated dependencies of the Purcell factor 
Fp (squares) and the radiative efficiency rj (diamonds) on the dipole 
position for a z-polarized dipole. The x-coordinate of the dipole is 
changed from -125 nm to 100 nm (measured from the position of 
the edge of the feed element), while the y- and z-coordinate are fixed. 
The vertical dashed lines mark the positions exactly of the edges of 
the director, feed, and reflector elements. The geometry is explained 
in the inset: the z-polarized dipole is shown as a black dot with arrow; 
it lies 10 nm above the feed's upper horizontal surface, is centered 
with respect to the feed element along v-direction, and shifted in x- 
direction, along the dotted horizontal line. 
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